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Abstract 

It this paper we present the results obtained with the newly created Standard SANC 
modules for calculation of the NLO QCD corrections to single top production processes 
in s and t channels at the partonic level, as well as top-decays. The main aim of these 
results is to prove the correct work of modules. A comprehensive comparison with results 
of the CompHEP system is given, where possible. These modules are intended to be used 
in Monte Carlo generators for single top production processes at the LHC. As in our 
recent paper, devoted to the electroweak corrections to these processes, we study the 
regularization of the top-legs associated infrared divergences with aid of the complex 
mass of the top quark. A comparison of QCD corrections with those computed by the 
conventional method is presented both for top production and decays. For s channel 
production we give an analytic proof of equivalence of the two methods in the limit of 
low top width. 

PACS: 12.15.Lk, 13.40.Ks, 14.65.Ha, 12.38.-t. 



1 Introduction 

There is continued interest in precision calculations of single top quark production cross 
sections for the Tevatron and LHC (see e.g. [I] and references therein). This is motivated 
to a large extent by the fact that this process is the only way of measuring the CKM 
matrix element |Vtb| directly, providing a sensitive test of the 3- generation scheme of the 
Standard Model, and that it is an important background to various processes beyond 
the SM, including Higgs boson production. 

Moreover, such events are already observed at the Tevatron experiments DO [21 13] 
and CDF [H [S] and at the LHC [B] . There is a need to prepare software for precision 
analysis of the high statistics samples of single top quark events in future measurements 
at the LHC, running at 7 and 14 TeV. 

Most of the theoretical work on single top production has been concerned with NLO 
QCD corrections (see e.g. [7], [8]), leading to the development of Monte Carlo generators, 
such as ZTOP [5], MC@NLO [TO] or SingleTop [H], incorporated in the standard LHC 
tools. 



iThis work is partly supported by RFFI grant Ar°10-02-01030-a 



In our previous paper [T2] we have presented the calculation of the cross sections with 
one-loop electroweak corrections, regulating the infrared divergences (IRD) associated 
with the top-quark leg by introducing the width of the top quark. 

In the present paper we extend that approach to the calculation of one-loop QCD 
corrections in the spirit of the complex- mass scheme f l3| . To our knowledge there is no 
other work utilizing this approach to the IRD regularization among the many papers 
dealing with higher order QCD corrections to single top quark production, see e.g. [H] 
and references therein, as well as papers devoted to the complex-mass scheme itself, 

e.g. [IS], dg. 

This paper should be read together with [12] , whose structure we have taken over here. 
In section [2] we recall the covariant amplitude for all processes thud — )■ and ibud —>■ 
and show the conversions to the processes of t and t decay and to s and t channel single 
top quark production, and point out the essential differences between the EW and QCD 
formulations of these processes. In section [3] we discuss in detail the regularization of 
single t and i quark production in the conventional approach — with zero top quark 
width — and in the complex top quark mass approach. We present explicit formulas 
derived for both cases. In section 2] we present the results of numerical calculations 
which show the stability against variations of the soft-hard separation parameter and 
of the top quark width. Finally, in section 5 we present our conclusions and outlook 
to further work on single top quark production within the framework of SANC, aimed 
at the creation of an MC generator that simultaneously takes account of NLO EW and 
QCD corrections in the so-called 5F-scheme [5]. 
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2 Covariant Amplitude 



In this section we proceed in the spirit of the "muhi-channel approach" developed in 
Refs. [IT] and [H], and fohow closely the presentation of our previous paper ^TT,. First, 
we consider annihilation into vacuum with all particles incoming. 

2.1 All particles incoming 

In QCD, there is no contribution from box diagrams such as those of [12] Figs. 1,2 which 
were appropriate for the EW case. Therefore in the present case we have a sum of two 
one- loop vertices instead of the two vacuum diagrams of [T2] : 

b d t U 





Figure 1: The tbud — >■ and tbud — >■ processes 

The Covariant Amplitudes (CA) of Figs. [T] are characterized by four different struc- 
tures and scalar Form Factors (FF) if the masses of the light quarks but not the b quark 
mass are neglected. A common expression for this CA in terms of four scalar form fac- 
tors, J^LL.RL.LD.RniSjt), was prescntcd in Ref. [ITj. We recall it here to introduce the 
notation: 



A — ie 



2 dvir{s) 



7m7+ 7p7+-7^tt(s, t) + 7^7_ (g) 7^7+ J'j^xls, t) 



+J+ (g) "ff,^+{-iDf,)J'LD{s,t) +7- (g) "ff,^+{-iDf,)J'Ro{s,t) 

= {pi -P2)m ■ 



where 

7± = (1±75) 
4-momentum conservation reads 



Pi + J32 + P3 + P4 = , 

and the invariants are defined by 
with 



s + i + u = rn^ -|- ml . 



Furthermore, 



2sf„ 



1 



(1) 
(2) 
(3) 

(4) 
(5) 
(6) 



where V = VtbVud is the relevant product of CKM matrix elements. The scalar form 
factors F are labeled according to their structures, see |17) . 
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2.2 Conversion to top decay 

The CA for the decay t(p2) b{pi) + it(p3) + d{p4) is derived from the left Fig. [T]by the 
following 4-monientuni replacement: 

Pi -Pi, P3 -P3, 

P2 P2, P4 -Pa- 

For the decay t(pi) — )■ h{p2) + dipz) + u{p4) it is more transparent to make the 
replacement from the right Fig. [T] 

Pi Pi, P3 -P3, 

P2 -P2, Pi -~P4- 

The corresponding two physical subprocess diagrams are those shown in Fig. 3,4 of [T^ . 



2.3 Conversion to s channel 

The CA for the s channel single-top production processes u{pi) + d{p2) ~^ ^(ps) + i{p4) 
is obtained from the left Fig. [1] by the replacement: 

Pi -P3, P3 Pi, 

P2 -Pa, Pa P2, 

and for the processes d(jii)+u{p2) — > t{p3) + b{p4) from the right Fig.[T]by the conversion: 



Pi -P3, P3 Pi, 

P2 -Pa, Pa -> P2- 

As a result we get the two physical subprocess diagrams shown in Fig. 5,6 of [T^ . 



2.4 Conversion to t channel processes 

In the CAs for the t channel single-top production processes b(pi) + u{p2) — >■ d(p3)+t{p4) 
and b{pi) + d{p2) — u{p3) + i(p4) it is convenient to make the replacement "in pairs". 
From the left Fig. [l]one may perform two 4- momentum replacements: 



Pi Pi, 

P2 -PA, 

P3 P2, 

Pa -> -P3, 



Pi Pi, 

P2 -PA, 

P3 -P3, 

PA P2, 



which give rise to two different physical t channel processes as shown in Fig. 7 of \TT . 

For the processes b{pi) + u{p2) t{p4) + d{p^) and b{pi) + d{p2) t{pi) + u(p3), 
the pair of replacements from the right Fig. (TJ 



Pi -Pa, 

P2 Pi, 

P3 -> -P3, 

PA P2, 



Pi ^ -PA, 

P2 Pi, 

P3 -> P2, 

PA -P3, 



gives the corresponding pair of symbolic diagrams for the two remaining t channel pro- 
cesses. They are shown in Fig. 8 of |12) . 
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3 Process ud tb at NLO QCD 

3.1 Notation and terminology 

After conversion of the vacuum annihilation diagrams of Fig. [1] to the s channel, we 
will have two physical annihilation diagrams for both t and t production (cf. symbolic 
diagrams in Fig. 5-6 of Ref. [E]). For our purpose it is sufficient to consider the top 
production NLO diagrams of Fig. 6 which, in turn, consist of ISR and FSR NLO con- 
tributions like the two fat vertices in diagrams in Fig. [l] standing for emission (and 
reabsorption) of a real (or virtual) gluon (cf. diagrams in Fig. 9,10 and 12 of Ref. [121). 

Since we are interested in the effects of the top quark width, the subject of our study 
will be only the symbolic FSR NLO diagram of Fig. [51 In the following, the label FSR 
will be assumed and omitted. It should be emphasized also that all the formulae of this 
Section refer to the limit rrib 0. 




Figure 2: The FSR in the du tb{g) processes 

We will use the following common terminology: 

— LO (or Born); 

— NLO=Virtual + Real; 

— Real=Soft + Hard; 

— SV =Soft + Virtual (with 2— >2 phase space); 

— NLO=SV+IIA (Hard Bremsstrahlung with 2—^3 phase space). 

3.2 Conventional approach 

In the conventional approach (stable particles) one introduces a soft-hard separation 
parameter 

UJ-l -^W, (7) 

with invariants s = ^ and s' — M"^ ^ for the radiative processes a + h ^ c + d + g. 
And for the total NLO cross sections we have 

The infrared divergence cancels inside cr^^ = 0"^+^ no matter how it is regularized 
(dimensionally or by a fictitious photon mass A), while the Co divergence cancels in the 

sum O-NLO ^ o-SV+HA^ go tj^^t ^NLO jg gj^j^g^ 

At the end of this paragraph, we list the individual contributions. 
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Soft+Virtual: 



'3mf 5 



-2 + (L, - 2) In (w) + ( - - Lti ) Lt 



2„ 2 S""n?^"i-'""r'"^"<^ 



+2U2{l)--Kmf-^{Ln+Lt) 
4 s 



(9) 



Hard: 



9 m^^ 



St 



{L,~2)\n{u) 



Lti L. 



(10) 



Born or LO : 

Here and below: 
asCf 



2 ^ 2Lti I Lti 



- - Lti ) Lt - 2Li2 (i; 



ii^i(2,s + m?). 



-i^ (st + 2TO^Lt) 



(11) 



fco = 



K = 



TT 



L, = 2 In 



St 



mtrrib 



s 



it =ln 



Lh = In 



Lti = In 1 - 



St — s — . 



(12) 



In the NLO cross section all divergences, including b mass collinear divergences, can- 
cel. The final expression is finite and very compact: 



NLO/ \ 7 ) Born/ \ 

a [s) = ko<(T (s) 



^-(Lt + l)Lti + 2Li2 

4 \ s 



1 j^rn'^ 
4 s2 



(5st (s + ™?) + 4:mfs) Lt - 3sf {Lti + 1) - 4m^st 



(13) 



3.3 Infrared regularization by the complex top quark mass 

The main aim of this paragraph is to show that the infrared regularization by the com- 
plex top quark mass leads exactly to the same final result, Eq. (jl3p . though the partial 
expressions for SV and HA have a completely different form. For the introduction to 
the problem, wc refer the reader to Section 3 of Ref. [IT. Here we again only list the 
individual contributions to the total cross sections. 
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Soft+Virtual: 



-1-(1-L,)ln( — 1 -ln(a>) 
rrit 



3 1 
- - 3Lii - -Lt I Lt 



-Lia 



3Li2(i; 



--Km^^^{Lti+Lt) 
4 



2Lti I Lti 

(14) 



Hard: 



^_!:!l + (l_£,)ln(Ii)+ln(c:,) 



St 



St 



mt 



"^Lti Lti 



2iti — I ^* 



- 3Li2 (1) 



(15) 



We note the presence of two regularization parameters: w is due to conventional 
treatment of infrared divergences associated with b quark legs, and \n(Tt/mt) is a new 
regularization parameter associated with the t quark leg, whose infrared divergence is 
regularized by its width. The individual (SV and HA contributions) have completely 
different form, but they sum up to exactly the same finite NLO expression (|13p in which 
all three types of divergent terms {\n{Tt/mt), In(aj), ln(mf,)) have canceled. 



3.4 Relevant AV functions 

The virtual contributions (due to vertex. Fig. 10, and the counter-term. Fig. 9, of 
Ref. [12) were computed using the standard Passarino-Veltman reduction [19j. The 
relevant Cq function 

Co {-rrit, -ml, Q^; rht, 0, rrib) , (16) 

with complex argument 



™2 



mf + At 



-imtTt, 



(17) 



is given by Eqs. (22)-(24) of 12^ and its limit at to;, — >■ by Eqs. (25)~(26) of that paper. 
Below we present its real part double limit: rub — >■ and Fj — >■ 0, which was used in the 
derivation of Eq. ([Tl)) : 



SRCo (-rrit, -nil, ~s; rht, 0, rrib) = — 

St 



mt 



Lti — Lt 



{Lb + Lt + 2Lti) (18) 



+LtiLt + -Lf - Li2 ( 1 - 



(l-^)+4Li2(l)|, 
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In the calculation of the counter-term, one meets the real part of the derivative of a Bq 
function: 

^[B'o{~mlO,mt)]^^\n(^] . (19) 
3.5 Differential in s' hard gluon radiation 

The hard gluon contributions, Eq. (fTU|) and Eq. (fT5|) . were derived by 4- fold integration in 
the 2—^3 process phase space over three angular variables and over the invariant s', which 
varies within the limits < s' < sj^^^^. In this paragraph we present the s'-integrands 
for the two approaches under consideration. 

3.5.1 Conventional approach 

In the conventional approach, the expression for (t^^(s, s') is well known: 

s') = fco|a^°"(s)(2 - L,,) (^i - - k'—^ (^^ - fciL.,) |, (20) 

where 

„' _ „' ™2 



s =s'-m^ (21) 

\mtmb J 

The w divergent term l/(s — s') after integration over s' gives rise to In (w) terms in 
Eq. (|10[) which cancel the corresponding terms in Eq. 

3.5.2 Infrared regularization by the complex top quark mass 

The expression for (t^^(s, s') is more complicated in this case: 



HA^ /\ 1 ) Born/ \ 

a (s,s ) = ko<(T (s) 



1 1 k+ , 1 

J + - -Ls' 

St s — s' 2 St 



;{h{s') + h{s'))--^h{s') 



-k'-^ [■^-k.L,. ] ). (22) 



Here, 

L' = - In 



s 



s — s' + l^t s — s' — l^t 

-ft = Ttmt . (23) 



The three non-trivial objects to be integrated are: 

h{s') = k+^{J^{Q,mt,V^,iit)). 
I^{s') = zfc-3(jr(0,mt,y?,z7t)). 



/3(,s') - l5(J,^(0,mt,Vs',*7t)), (24) 
7t 
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where 



/ ijts' + (s ~ s')m^ 
\ ijts' + (s — s')s' 



and "0" in the argument hst stands for the hmit mi, — > 0. 
Taking the integrals 



his)= f his') ds', ^ = l,2,3, 



(25) 



(26) 



we get 



his) = In — 



ln(lL)-21n(-f^ 



his) 



1 



-St In 



nit 



St In 



'St 



In I 



(27) 



Substituting Eqs. (H?]) into Eq. dm we arrive at Eq. ([TS 
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4 Numerical Results 



In this section we present the SANC results for the cross sections of the single top quark 
production processes and for the top quark decays. The tree level contributions, both 
Born and single real gluon emission, are compared with CompHEP versions v. 4. 4. 3 pD] 
and v.4.5.1 [21 . All numerical results for this section were produced with the standard 
SANC INPUT working in the a{0) EW scheme. 
The SANC input parameters set: 



Gp 


= 1.16637 X 10-5 GeV"^ 






a(0) 


= 1/137.035999, 




= 0.107 


M„ 


= 80.403 GeV, 




^ 2.141 GeV, 




= 91.1876 GeV, 




= 2.4952 GeV, 


M„ 


= 120 GeV, 


r* 


= 1.5517 GeV, 


niu 


= 62 MeV, 


rrid 


= 83 MeV, 


rric 


= 1.5 GeV, 


rus 


= 215 MeV, 


mt 


= 174.2 GeV 


rUb 


= 4.7 GeV, 


\Vud\ 


= 1, 


\Vcs\ 


= 1, 


\Vus\ 


- 0, 


\Vcd\ 


= 0, 


\Vtb\ 


= 1. 







Numbers used for the comparison with CompHEP for the decay channels were produced 
with CompHEP version 4.4.3 and with version v.4.5.1 for the production channels. For 
this comparison we use the standard CompHEP setup and the a(0) EW scheme, however, 
we used as = 0.12201 {Q = Mz) and mu = md = 66 MeV. Moreover, in both codes we 
use the "fixed t width scheme" . 

The structure, notation and terminology used of this Section are very similar to those 
of the corresponding Section 4 of Ref . [12] . We refer the reader to the paragraphs that 
follow Eq. (33) and do not repeat here Eqs. (34)-(35) and the surrounding text of the 
latter paper. 

4.1 Decay channels t ^ b + u + d 

4.1.1 SANC CompHEP comparison 

In Table [T] we present results of the SANC-CompHEP comparison of radiative decay width 
phaid MeV for two cuts on the gluon energy: Eg > 5, 10 GeV (in the top rest 
frame) and for two options: Tt — 0{y^ 0) and six channels pid=19-24 for the case of 
SANC (pid=19,21,23 for i decays, pid=20,22,24 for t decays). 

The main aim of Table [T] is to demonstrate, for the first time, that there is good 
agreement between SANC and CompHEP not only for Ff — 0, but also for the Ff 7^ 
options (Ft 7^ always refers to the value in Eq. (|28p '). 

4.1.2 An internal SANC comparison 

In SANC there are two modules for the computation of hard gluon bremsstrahlung pro- 
cesses: fully differential (5d), to be used in Monte Carlo simulations, and two dimensional 
(2d) in invariant variables s = —{pt — PbY ^'^'^ ^' — ~{Pu + PdY , being analytically in- 
tegrated over all angular variables. In order to check numerically the correctness of 
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Eg, GeV 


5 


10 




r* = 


^0 


-0 


r* ^0 


t^b + e+ + i?e+g 


CompHEP 


72.29(1) 


71.80(2) 


47.24(2) 


47.13(2) 


SANC , pid=19 


72.28(1) 


71.79(1) 


47.20(1) 


47.10(1) 


SANC , pid=20 


72.28(1) 


71.79(1) 


47.20(1) 


47.10(1) 


t ^ b + fi+ + + g 


CompHEP 


72.29(2) 


71.80(2) 


47.24(2) 


47.13(2) 


SANC , pid=21 


72.28(1) 


71.79(1) 


47.20(1) 


47.10(1) 


SANC , pid=22 


72.28(1) 


71.79(1) 


47.20(1) 


47.10(1) 


t^b+u+d+g 


CompHEP 


1296.9(6) 


1295.6(5) 


820.3(3) 


819.7(3) 


SANC , pid=23 


1296.3(1) 


1294.9(1) 


819.8(1) 


819.5(1) 


SANC , pid=24 


1296.3(1) 


1294.9(1) 


819.8(1) 


819.5(1) 



Table 1: SANC-CompHEP comparison of decay width r in MeV in the a scheme for 
two options: Ft = 0(7^ 0); and for six channels pid=19-24. 

the analytic angular integration, we present Table [5] which shows the self-consistency of 
SANC calculations of the radiative contribution. Note that here Eg is the cut in the s' 
compound rest frame: this is why the numbers in Tables [T] and [5] are different. 



Eg, GcV 


5 


10 




Ft =0 




r* = 


Tt^O 


t-^b + 


/i+ + i^^ + .g 


SANC-5d, pid=22 


96.446(1) 


94.459(1) 


67.492(1) 


67.085(1) 


SANC-2s', pid=22 


96.446(1) 


94.459(1) 


67.493(1) 


67.085(1) 


t ^ 5^ 


- u + d + g 


SANC-5d, pid=24 


1245.6(1) 


1239.6(1) 


763.09(1) 


761.87(1) 


SANC-2s', pid=24 


1245.6(1) 


1239.6(1) 


763.11(1) 


761.88(1) 



Table 2: SANC comparison of decay width r^'*''^ in MeV in the a(0) scheme for two 
options: Ft = 0(7^ 0); two variants: 5d and 2s'; and for two channels pid=22 and 24. 

4.1.3 One-loop decay width 

The numbers of this Section are produced with SANC setup (|28l) in the Q!(0) EW param- 
eterization scheme. We present results for two processes: pid=22 and 24. 

• pid — 22. Process t — > 6 + /i+ + i/^ 

In Table [3] we show the lowest-order width r^°''" as a function of the 6 quark mass. 
As is seen, the effect of a finite b quark mass is tiny. This justifies considering the limit 
TOft — 0, in which the formulas are much shorter and which is convenient to study the 
validity of the KLN theorem in NLO approximation. 
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nib, GeV 


4.7 


1.0 


0.1 


r^°™, GeV 


0.149094(1) 


0.149466(1) 


0.149483(1) 



Table 3: The total lowest-order width r " in GeV as function of mf,. 

Table m demonstrates the stability of one- loop corrected quantities to the variation of 
the soft-hard separation parameter Eg. 



Tt=0 


Eg, GeV 


10-' 


10-^ 


10-^ 


pi — loop 


0.13651(1) 


0.13652(1) 


0.13651(1) 


d% 


-8.440(1) 


-8.437(2) 


-8.439(9) 




pi — loop 


0.13735(1) 


0.13734(1) 


0.13735(1) 


d% 


-7.880(2) 


-7.883(2) 


-7.877(8) 



Table 4: The total one-loop corrected width '°°p in GeV and relative one-loop 
correction S = (ri-i°°P - rB°™)/ri-i°°P for Eg = lO^i, IQ-^, lO'^ GeV and TOfc=4.7 
GeV. 

As is seen from Table SI there is good Eg stability for Eg < 10^^ GeV. However, 
there is a rather big difference of the results for the one-loop corrected width r^~'°°P and 
the relative QCD correction S of the partial decay width between the calculation with 
the = and the Fi 7^ options, where in the latter case we use the fixedi scheme 
defined in Ref. 12 . This difference gets smaller with formally decreasing Ff as Table [S] 
illustrates. 

Table [S] demonstrates the two dimensional convergence of r^^'°°P and S when both 
mb and Tt go to zero. Note that with decreasing Fj the convergence in mi, improves. 



pi-ioop^ GeV 


rt/io'\ N 





1 


2 


Ff = 


ruh, GeV 










4.7 


0.137345(2) 


0.136594(2) 


0.136518(3) 


0.136515(3) 


1.0 


0.137919(2) 


0.136840(2) 


0.136726(3) 


0.136717(3) 


0.1 


0.138479(2) 


0.136900(2) 


0.136738(3) 


0.136723(5) 


5,% 


Tt/W\ N 





1 


2 


Ff = 


nib, GeV 










4.7 


-7.881(1) 


-8.384(1) 


-8.436(2) 


-8.437(2) 


1.0 


-7.726(1) 


-8.448(1) 


-8.523(2) 


-8.530(2) 


0.1 


-7.362(1) 


-8.418(1) 


-8.526(2) 


-8.537(3) 



Table 5: The total one-loop corrected width r^^'°°P in GeV and corresponding relative 
one-loop correction 6 in % as function of the b quark mass and of Ff with soft-hard 
separation parameter Eg = 10^^ GeV. 
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• pid = 24. Process t ^ b + u + d 

Here we present the corresponding set of Tables for the quark decay mode of the top 
quark. 

The total lowest-order width (Table [S]) also shows weak sensitivity to variations of 
the b quark mass. 



nib, GeV 


4.7 


1.0 


0.1 


r^°™, GeV 


0.447284(1) 


0.448399(1) 


0.448451(1) 



Table 6: The total lowest-order width r " in GeV as function of nib- 
The stability against variation of Eg is shown in Table [71 





Eg, GcV 


10-' 


10-^ 


10-^ 


pi— loop 


0.3594(1) 


0.3581(1) 


0.3584(4) 


S % 


-19.65(1) 


-19.94(1) 


-19.89(6) 


7^ 0, fixedi 


pi— loop 


0.3619(1) 


0.3606(1) 


0.3607(3) 


6% 


-19.10(1) 


-19.37(1) 


-19.35(7) 



Table 7: The total one- loop corrected width '°°p in GeV and relative one-loop 
correction 5 in % for Eg ^ 10-\ 10"^ lO'^ GeV and mb^A.l GeV.. 

Similar conclusions as drawn after Table S] are valid in this case. Finally, Table [51 
showing the two dimensional limit, demonstrates similar behavior as Table [SI for the 
semi-leptonic top decay. 



pi-ioop^ GeV 


Ft/10^\ N 





1 


2 


Ft =0 


nib, GeV 










4.7 


0.36054(8) 


0.35830(8) 


0.35808(8) 


0.35812(6) 


1.0 


0.36210(8) 


0.35890(8) 


0.35863(8) 


0.35870(6) 


0.1 


0.36380(8) 


0.35898(8) 


0.35854(8) 


0.35861(7) 


s,% 







1 


2 


Ft =0 


TTib, GeV 










4.7 


-19.39(2) 


-19.89(2) 


-19.94(2) 


-19.94(1) 


1.0 


-19.25(2) 


-19.96(2) 


-20.02(2) 


-20.00(1) 


0.1 


-18.88(2) 


-19.95(2) 


-20.05(2) 


-20.03(2) 



Table 8: The total one-loop corrected width F^ '°°p in GeV and corresponding relative 
one-loop correction 5 in % for Eg = 10~^ GeV as function of the b quark mass and of Ft. 
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4.2 s channel 

Here we consider s channel processes: 

• pid=25, J + u — > 6 + t and pid=26, u + d — > t + 5 
but it is sufBcient to consider one channel, say pid=25. 

4.2.1 SANC-CompHEP comparison 



Vs, GeV 


200 


1000 


7000 


-0 


CompHEP 


0.72930(2) 


1.6340(1) 


0.071223(1) 


SANC-4d 


0.72927(1) 


1.6341(1) 


0.071229(1) 


SANC-s' 


0.72928(1) 


1.6338(1) 


0.071227(1) 




CompHEP 


0.72638(2) 


1.6323(1) 


0.070946(1) 


SANC-4d 


0.72636(1) 


1.6322(1) 


0.070942(1) 


SANC-s' 


0.72636(1) 


1.6322(1) 


0.070942(1) 



Table 9: Comparison of the cross section cr^"'^(\/l, Fj) in pb for three cms energies 
and two width options Ft — 0). for pid = 25, i.e. for process d + u b + t; the Eg 
cut is equal to 2 GeV. 

As is seen from Table [3J there is very good agreement of numbers obtained from 
SANC and CompHEP within the statistical errors for two considered width options and 
three cms energies. In this case, the two widths options agree well for all three energies. 

4.2.2 One-loop corrections 

The analogue of Tables 4-9 of Ref. [12] is now represented by one joint Table [TOl showing 
the stability of one-loop corrected QCD cross sections (t^~'°°p and relative QCD RC 5 
against variation of the soft-hard separation parameter, tj, and the difference between 
the two options: Ft = 0(^ 0). 

As is seen, there is good ui stability in the interval ui — 10^^ -f- 10^^; also, the 
subtracted quantities and for all three cms energies and for both width options, 
r* — 0(7^ 0), are seen to be independent of the light quark masses. 

One can see that the difference between the two width options is of order of few per 
mille in absolute deviation, rapidly decreasing with increasing cms energy. We therefore 
conclude that one may use the usual infrared regularization for s channel processes as 
was the case for EW corrections, see text after Table 6 of Ref. |12j . 

In Table [TT] we demonstrate two dimensional convergence of tr^^ and when both 
mb and Ft go to zero. Note that, similarly to top decay cases, with lowering of Ft the 
convergence in nih improves. 

The agreement of af^^ and in the double limit (numbers for Tt/lO* and Ft = 0) 
illustrates our analytic considerations described in Section [3] 
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yi, GeV 


200 


1000 


7000 




308.09055(1) 


105.97718(1) 


2.2352950(1) 








10"'^ 


lO"** 


Tt = 0, ruq 


Ijl-loop 


-374.05(1) 


-374.07(1) 


0.36446(1) 


364.46(1) 


17.599(1) 


17.599(1) 


5,% 


-221.41(1) 


-221.42(1) 


243.90(1) 


243.90(1) 


687.34(1) 


687.34(1) 




525.34(1) 


525.37(1) 


0.17057(1) 


170.57(1) 


8.2097(1) 


8.2097(1) 




70.51(1) 


70.53(1) 


60.95(1) 


60.95(1) 


267.28(1) 


267.28(1) 


Ft =0,m,/10 




525.34(1) 


525.37(1) 


0.17057(1) 


170.57(1) 


8.2097(1) 


8.2097(1) 




70.51(1) 


70.52(1) 


60.95(1) 


60.95(1) 


267.27(1) 


267.28(1) 






-372.85(1) 


-372.87(1) 


0.36455(1) 


364.55(1) 


17.601(1) 


17.601(1) 


^, % 


-221.02(1) 


-221.03(1) 


243.99(1) 


243.99(1) 


687.40(1) 


687.40(1) 




526.54(1) 


526.57(1) 


0.17066(1) 


170.66(1) 


8.2111(1) 


8.2111(1) 




70.90(1) 


70.91(1) 


61.03(1) 


61.03(1) 


267.34(1) 


267.34(1) 


Ft ^0,mg/10 




526.54(1) 


526.57(1) 


0.17066(1) 


170.66(1) 


8.2111(1) 


8.2111(1) 




70.90(1) 


70.91(1) 


61.03(1) 


61.03(1) 


267.34(1) 


267.34(1) 



Table 10: The total lowest-order cross section ct^""^", one- loop cross section '°°p and 
total one-loop corrected MS subtracted quantities ij'^^ and corresponding 5 and in 
% at w = 10-^ 10-6 for ^ = 200, 1000, 7000 GeV. 



af^, pb 


Ft/10^\ N 





2 


4 


Ft =0 


mb, GeV 










4.7 


0.52654(1) 


0.52530(1) 


0.52534(1) 


0.52534(1) 


1.0 


0.53099(1) 


0.52783(1) 


0.52776(1) 


0.52778(1) 


0.1 


0.53386(1) 


0.52791(1) 


0.52784(1) 


0.52784(1) 


0.01 


0.53656(1) 


0.52793(1) 


0.52783(1) 


0.52785(1) 




Tt/W\ N 





2 


4 


Ft =0 


nib, GeV 










4.7 


70.90(1) 


70.50(1) 


70.52(1) 


70.51(1) 


1.0 


69.65(1) 


68.64(1) 


68.62(1) 


68.62(1) 


0.1 


70.44(1) 


68.54(1) 


68.52(1) 


68.52(1) 


0.01 


71.. 30(1) 


08.33(1) 


G8.32(l) 


08.32(1) 



Table 11: The total one- loop corrected MS subtracted cross sections pb and 

corresponding relative one-loop correction in % for \/I=200 GeV (at the parton 
level) and for w = 10"^ as function of the b quark mass and of Tf 
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4.3 t channel 



4.3.1 SANC CompHEP comparison 

This is the most complicated case: t channel cross sections show up bad statistical 
convergence. For this comparison we use CompHEP v. 4. 5.1 and its setup, but with non- 
zero masses of the u and d quarks (accessed via be — ^ tsj channel). For the Tables of 
this subsection we used rrig = rriu — ma = 66 MeV and lOmg, i.e. m„ = to^ = 660 MeV. 
The cut on the cms gluon energy was Eg > 2 GeV, and as = 0.12201 {Q — Mz). As 
in [12], rows marked "SANC(S)" ("S"hort massive case) were computed retaining niq or 
lOwg only in fermion propagators radiating a gluon, while "SANC(F)" ("F"ully massive 
case) means that light quark masses were kept everywhere. 

• pid = 27: Process 6 + u — > rf + t 



Vs/GeV 


200 


1000 


7000 


Vt = 0, niq 


CompHEP 


17.908(1) 


415.04(6) 


617.40(6) 


SANC(S) 


17.904(1) 


416.39(1) 


618.94(1) 


Vt ^ 0, niq 


CompHEP 


17.761(1) 


404.29(6) 


574.02(6) 


SANC(S) 


17.759(1) 


405.62(1) 


575.61(1) 


SANC(F) 


17.760(1) 


405.63(1) 


575.62(1) 


Vt ^ 0, IQrUq 


CompHEP 


11.792(1) 


286.79(7) 


400.29(2) 


SANC(F) 


11.789(1) 


287.08(1) 


400.48(1) 



Table 12: Comparison of the cross section cr^^^'^{-/s) in pb for the process h + u ^ t + d 
for three cms energies; three options of four: {Vt — 0(^ 0)) (8) {ruq, lOmq). 

• pid — 28: Process 6 + d — J- u + t 



VI/GeV 200 


1000 


7000 


Vt = 0, TOg 


CompHEP 


10.335(1) 


388.74(7) 


615.5(1) 


SANC(S) 


10.332(1) 


389.87(1) 


616.86(1) 


Vt ^ 0, niq 


CompHEP 


10.245(1) 


378.21(7) 


572.0(1) 


SANC(S) 


10.241(1) 


379.34(1) 


573.58(1) 


SANC(F) 


10.241(1) 


379.35(1) 


573.59(1) 


Vt 0, IQmq 


CompHEP 


6.641(1) 


266.75(2) 


398.68(4) 


SANC(F) 


6.640(1) 


266.83(1) 


398.83(1) 



Table 13: The same as Table [T^l but for the process h + d ^ u + t. 
The results of Tables [T^HT51 are qualitatively the same and may be discussed together. 
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There is satisfactory agreement for all options only near the reaction threshold. At 
higher energies SANC/CompHEP agree poorly for the niq option and much better for the 
lOrUq option. This allows us to draw conclusions similar to those which were drawn at 
the discussion of the results of analogous Tables 10-11 of Ref. |12j . 

4.3.2 One-loop QCD corrections 

The numerical results for this subsection were again produced with the aim to demon- 
strate the stability of one-loop corrected QCD cross sections o-1~'°°p and relative QCD 
RC S against variation of the soft-hard separator uj and to study the difference between 
the two options: Ft = 0(7^ 0). 

• pid = 27. Process 6 + u — > rf + t 

Table [His the counterpart for QCD of the six Tables 12-17 of Ref. HI]. 



y/s, GeV 


200 


1000 


7000 


cr^°™, pb 


7.3551153(1) 


48.993407(2) 


50.824228(2) 




10-'^ 


10-^ 


10-^ 


10-^ 


10-^ 


io-« 


Ft = 0, ruq 


0-l-loop 


-5.741(1) 


-5.742(1) 


43.67(1) 


43.66(1) 


50.57(2) 


50.57(2) 


s,% 


-178.06(1) 


-178.07(1) 


-10.90(2) 


-10.88(3) 


-0.50(3) 


-0.48(3) 




10.215(1) 


10.215(1) 


47.67(1) 


47.67(1) 


50.71(2) 


50.72(2) 




38.88(1) 


38.88(1) 


-2.70(2) 


-2.69(3) 


-0.23(3) 


-0.20(4) 


Ft =0,m,/10 




10.212(1) 


10.210(1) 


47.66(2) 


47.66(2) 


50.71(2) 


50.72(3) 




38.88(2) 


38.86(2) 


-2.73(4) 


-2.72(4) 


-0.22(5) 


-0.21(6) 




(jl-loop 


-5.573(1) 


-5.574(1) 


43.94(1) 


43.95(1) 


50.87(2) 


50.86(2) 


s, % 


-175.77(1) 


-175.78(1) 


-10.29(3) 


-10.28(2) 


0.097(29) 


0.063(34) 



Table 14: The total lowest-order cross section cr^"""", one- loop cross section '°°p and 
total one-loop corrected MS subtracted quantities (all in pb) and corresponding S 
and in % for w = 10-^ 10"^ and = 200, 1000, 7000 GeV. 

Stability in w is seen to hold in all considered variants; niq independence is clearly 
seen for the Fj = option. However, contrary to EW corrections, the dependence on 
Ft itself for non-subtracted quantities (compare two 6 rows) is much more pronounced, 
varying from several % near the threshold to half a percent at higher energies. Since it 
is not known how to realize subtraction for F^ ^ option, this difference may be treated 
as a theoretical uncertainty of NLO calculations which is, however, much smaller than 
the estimate of NNLO contributions [12] • 

For t channel processes, we show only one-dimensional convergence in the limit Ft — > 
for y/s = 200 GeV where our NLO correction is unphysical (< —100%), see Table [T5l 
(We recall, that in this case we do not have an analytic proof of convergence similar to 
the proof shown in Section [3] for the s channel). So, it should be considered as only a 
formal, numerical illustration of convergence. 
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= 200 GeV, cr"°™ = 7.3551153(1) pb 




^l-loop^ pb 


6,% 





-5.5730(6) 


-175.77(1) 


2 


-5.7392(6) 


-178.03(1) 


4 


-5.7419(7) 


-178.07(1) 


= 


-5.7411(6) 


-178.06(1) 



Table 15: The total lowest-order cross section (t^°'", the total one-loop corrected tr^ 
in pb and corresponding (5 in % at = 200 GeV and Tt/lO^, iV = 0, 2, 4. 



• pid = 28. Process b + d ^ u + t 

Table [m is the counterpart for QCD of the six Tables 18-22 of Ref. [T^. 



y/^, GeV 


200 


1000 


7000 


cr"°™, pb 


4.49579065(2) 


46.6869560(1) 


50.72449(1) 




10"^ 10-" 


10-" 10-" 


10-" 10-" 


Tt — 0, rriq 


^l-loop 


-3.2657(4) 


-3.2666(6) 


39.71(1) 


39.70(1) 


50.22(2) 


50.22(2) 


<5,% 


-172.64(1) 


-172.66(1) 


-14.95(2) 


-14.96(3) 


-0.99(3) 


-1.00(4) 




6.5587(4) 


6.5586(6) 


45.29(1) 


45.29(1) 


50.53(2) 


50.52(2) 




45.88(1) 


45.88(1) 


-2.98(2) 


-2.99(3) 


-0.38(3) 


-0.39(4) 


Ft =0,m,/10 




6.5564(7) 


6.5553(9) 


45.28(2) 


45.28(2) 


50.51(3) 


50.53(3) 




45.91(2) 


45.88(2) 


-3.01(4) 


-3.02(4) 


-0.42(5) 


-0.39(6) 


Tt ^ 0, ruq 


^l-loop 


-3.1526(4) 


-3.1532(5) 


39.99(1) 


39.99(1) 


50.52(1) 


50.51(2) 


<5, % 


-170.12(1) 


-170.14(1) 


-14.35(2) 


-14.34(3) 


-0.40(3) 


-0.41(3) 



Table 16: The same as Table [H but for pid = 28. Process b + d~>u + t. 



= 200 GeV, cr"°™ = 4.495790646(3) pb 




^1 — loop 


S 





-3.1526(4) 


-170.12(1) 


2 


-3.2644(4) 


-172.61(1) 


4 


-3.2661(5) 


-172.65(1) 


Ft -0 


-3.2657(4) 


-172.64(1) 



Table 17: The same as Table [13 but for pid = 28. Process b + d u + t. 

From Tables [T5HT71 the same conclusions may be drawn as from Tables [TH - fTSl We 
will not repeat them here. 
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5 Conclusions and Outlook 



In two papers, this one and Ref. |12) . we have described the implementation into the 
SANC framework of the complete one-loop QCD and EW calculations, including hard 
bremsstrahlung contributions, for the processes of top quark decays and of s and t channel 
production, the latter two at the partonic level. Within the SANC framework, we have 
created the standard FORM and FORTRAN modules, see, Ref. [23], compiled into a 
package sajic_cc_vl .40 which may be downloaded from the SANC project homepages 

The essentially new aspect of these two papers is the study of regularisation of the top 
quark-photon and -gluon infrared divergences with aid of the complex mass of the top 
quark. A comparison of these NLO corrections computed within this approach with those 
computed by the conventional method showed a difference which is not negligible, though 
not exceeding ~ 1% for both types of corrections and for three considered processes. 
Although, the introduction of the non-zero top width within the fixed width scheme 
violates gauge-invariance, it is numerically not so important being of the order 0{Tt/mt); 
it restores in the limit Ft — )■ 0. 

The emphasis of these papers is to be assured of the correctness of our results. We 
observe the independence of the form factors on gauge parameters, checked the stability 
of the result against variation of the soft-hard separation parameter lj and the indepen- 
dence of the MS subtracted quantities off the initial quark masses which is crucial for 
calculations at the hadronic level. 

As has become SANC standard, we have compared our numerical results with other 
independent calculations wherever possible. For the decay channels it was done in our 
previous papers (both EW [33] and QCD [IS]), showing good agreement. As usual, the 
Born level and the hard gluon contributions of all three channels were checked against 
the CompHEP package. We found very good agreement for both Ft = 0(7^ 0) options. 

As far as the comparison of EWRC for the production processes is concerned, it was 
described in our previous paper Ref. [12j and we found no suitable results in the literature 
on QCD corrections at the partonic level. 

As a first step toward comparison of QCD calculations at the hadronic level we per- 
formed a triple comparison of the LO total cross sections of s and t channel single top 
production computed by SANC/CompHEP/MCFM with the same set of PDF. We found a 
satisfactory agreement between SANC and CompHEP numbers. The comprehensive com- 
parison of the results at the hadronic level is the subject of ongoing work and a further 
publication. 

All the calculations were done using a combination of analytic and Monte Carlo 
integration methods which will make it easy to impose experimental cuts in future calcu- 
lations for pp collisions at the hadronic level. The results presented in these two papers 
lay a solid base for future extensions of calculations for the single top production chan- 
nels at the LHC with subsequent decay in the cascade (factorized) approximation (see 
e.g. [57]) which simultaneously take account of NLO EW and QCD corrections. 
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